Real-space renormalization group flow in quantum impurity systems: 
local moment formation and the Kondo screening cloud 
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The existence of a lengtli-scale ~ (with Tk the Kondo temperature) has long been 

predicted in quantum impurity systems. At low temperatures T <C Tk, the standard interpretation 
is that a spin-i impurity is screened by a surrounding 'Kondo cloud' of spatial extent £^k- We argue 
that renormalization group (RG) flow between any two fixed points (FPs) results in a characteristic 
length-scale, observed in real-space as a crossover between physical behaviour typical of each FP. 
In the simplest example of the Anderson impurity model, three FPs arise; and we show that 'free 
orbital', 'local moment' and 'strong coupling' regions of space can be identified at zero temperature. 
These regions are separated by two crossover length-scales ^lm and ^x, with the latter diverging as 
the Kondo effect is destroyed on increasing temperature through Tk- One implication is that moment 
formation occurs inside the 'Kondo cloud', while the screening process itself occurs on flowing to the 
strong coupling FP at distances ~ ^k- Generic aspects of the real-space physics are exemplified by 
the two-channel Kondo model, where £^k now separates 'local moment' and 'overscreening' clouds. 

PACS numbers: 72.15.Qm, 73.63.Kv, 75.20.Hr 
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I. INTRODUCTION 

Most fundamental aspects of the Kondo effect are by 
now very well understood,^ with various detailed theoret- 
ical predictions having been confirmed directly by exper- 
iments on impurity systems or quantum dot devices. 
Key insights into the underlying physics have been pro- 
vided by the renormalization group (RG) concept, where 
progressive reduction of the temperature or energy scale 
results in RG flow between 'fixed points' (FPs) that can 
be easily identified for a given model. In particular, the 
low-temperature/energy physics is governed by RG flow 
between two fundamental FPs, with universality arising 
in terms of the crossover energy scale. 

The Anderson impurity model (AIM) captures many 
generic aspects of quantum impurity physics. Three FPs 
(and hence two crossover energy scales) arise: charge fluc- 
tuations on the impurity or quantum dot dominate at the 
'free orbital' (FO) FP, which describes the high-energy 
physics. But due to strong impurity electron correlations, 
a spin-i 'local moment' (LM) forms below T ^ Tlm (typ- 
ically a high-energy scale) . RG flow from this LM FP to 
the 'strong coupling' (SC) FP occurs on a much lower en- 
ergy scale Tk, and determines the universal low-energy 
behaviour. The SC FP itself describes the ground state in 
which the incipient moment associated with the LM FP 
is dynamically screened by conduction electrons, which 
together form a many-body 'Kondo singlet' state. ^ 

Surprisingly however, the underlying physics in real- 
space is still somewhat controversial. Various theoreti- 
cal studies^^^^ have examined real-space behaviour in a 
number of impurity models. While they do appear to 
have established the existence of a Kondo length-scale, 
= /ksTK (with vp the Fermi velocity), the na- 
ture of the two spatial regions separated by it remains 
unclear. Certain quantities have been studied in the uni- 
versal scaling regime, ^"^^ but the full and exact evolution 
of real-space quantities — and a satisfactory unifying in- 



terpretation of the results — has not yet been provided. 

The prevailing view^*^ is that at low temperatures 
T ^ Tk , a spin- i impurity forms a singlet by binding to 
a surrounding 'Kondo cloud' of spatial extent ^k- This 
has an intuitive appeal from the perspective of Fermi 
liquid theory, where one imagines that the impurity is 
'invisible' to conduction electrons outside the screening 
cloud. ^'^ But no clear evidence in support of the screening 
cloud scenario has so far emerged. Very recently, refine- 
ment of scanning tunneling spectroscopy techniques re- 
ported in Ref. 16 has shown for the first time that Kondo 
signatures in the local density of states (LDOS) appear 
away from an impurity. Experimental verification of the 
elusive Kondo length-scale — and resolution of the real- 
space debate — may now finally be within reach. 

Motivated in part by these developments, we consider 
here the wider question of how RG flow manifests in 
real-space. While necessarily perturbative, a simple in- 
tuitive picture is suggested by well-known scaling ar- 
guments applied directly to real-space, where notions 
of RG and universality again arise. One might antici- 
pate then that this universality could show up in certain 
real-space quantities. Indeed, some concrete evidence 



Figure 1: Three physical re- 
gions arise at T = in the An- 
derson impurity model, cor- 
responding to the FO, LM 
and SC FPs (a). The 
crossover length-scale £_k as- 
sociated with flow to the SC 
FP diverges as the Kondo ef- 
fect is destroyed on increasing 
temperature through Tk (b), 
while ^LM diverges for T 3> 
Tlm (c). The spin-i moment 
forming at ^lm is screened on 
flowing to the SC FP at . 
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to support this general expectation has been provided 
theoretically^"^^ if not experimentally. 

Encouraged by these results, we show in this paper that 
motion away from an impurity corresponds to RG flow. 
In fact, we demonstrate explicitly and exactly that signa- 
tures of the entire RG flow between all FPs (whether uni- 
versal or non-universal) must appear fully in real-space. 
Just as flow between a pair of FPs results in a character- 
istic energy scale, we show that one can similarly identify 
the corresponding length-scale, apparent in real-space as 
a crossover between physical behaviour typical of each 
FP. In particular, this implies that length-scales can be 
attributed both to moment formation and moment com- 
pensation. In real impurity or quantum dot systems, the 
moment screened by the Kondo effect at low tempera- 
tures does not reside on the impurity or dot itself, but is 
an extended object which develops in real-space due to 
electron correlations. 

Further insight is gained by considering the effect on 
the real-space physics of increasing temperature. RG flow 
is cut off at an energy scale ~ T; likewise RG flow is 
arrested on the thermal length-scale oc 1/T. Conse- 
quently, clouds associated with low-energy FPs reached 
at larger distances 'evaporate' successively on increas- 
ing temperature. Our refined real-space Kondo screening 
scenario is summarized schematically in Fig. 1. 



II. ANDERSON IMPURITY MODEL IN 
REAL-SPACE 

For simplicity and concreteness, we con- 
sider first a semi-infinite tight-binding chain, 

-ffhost = tJ2aJ2Zoi4,a(^^+l,cr +H.c), with a single 
Anderson impurity tunnel-coupled to one end: 



-^And = Hhost + <^nd + Ufld.^fld.X + V 



(4co,. + H.c), 



(1) 

where fid = J^a- ^d,,a = J2cr '^a-d.a- '^^ the total number 
operator for the impurity, and the hybridization strength 
is r = TTpV^ (with p = l/(7rt) the Fermi level density of 
states and 2t = 1 the half-bandwidth). As in previous 
theoretical work,^"^ we focus here on the 'excess' charge 
density due to the impurity, obtained via 
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dij f{Lj,T) AGrA^,T), (2) 



with f{uj,T) = [1 +exp(a;/T)] ^ the Fermi function and 
AGr,aiu},T) = Gr,a{uJ,T) - G^{uj) the difference in the 
site-r Green function with and without the impurity. 

The free conduction chain is initially at half-filling, 
and the corresponding no-impurity propagators G'^{uj) 
are simple objects, describing non-interacting electrons. 
By contrast, the behaviour of the full Green function 
Gr,a{'jJ,T) = {{c^ ^; cl ^))i^^T IS deeply non-trivial due 
to electron correlations on the impurity which are ul- 
timately responsible for local moment formation and 
thence the Kondo effect. But these correlations are lo- 
cal, and so the subtle real-space behaviour of Gr,a{oj,T) 
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Figure 2: Comparison of the spectrum vs inverse frequency 
(upper panel) and normalized excess charge density vs dis- 
tance (lower panel) for the Anderson model at T = 0. Plotted 
for common lO^F = 3 and lO^e = -0.8, varying lO^U = 1, 
1.125 and 1.25 (solid, dotted and dashed lines), chosen to give 
exaggerated energy- and length-scale separations for clarity. 



and An{r, T) can be determined from purely local impu- 
rity quantities. This is most clearly seen using equations 
of motion,^ whence for the simple geometry considered 
here, one readily obtains 



AGrA^,T) = [\XiLj,r)]^ T„{^,T). 



(3) 



All the real-space information about the lattice is con- 
tained in the simple function X{La,r) = [tG^{ijj)Y'^^ , 
where GI{lu) = (w - Vw^ - At'^)/2t^ is the free r = 
Green function. The spatial dependence of AGr,(T(w,r) 
thus enters only through the power of G^{ijj), reflecting 
physically the fact that electrons must hop {r-\-l) times to 
get to the impurity. T„{ijj, T) = V^Gd,a (w, T) is the usual 
scattering t-matrix, with Gd,cr{'^,T) = {{d^;d'l))^,T the 
impurity Green function. In the trivial non-interacting 
limit [U ~ 0), the impurity merely gives rise to po- 
tential scattering, and the familiar Friedel density os- 
cillations are expected.^* In this case G^J^^{ljj) = [uj — 
e — V'^Gq{uj)]^^ independent of temperature. Far from 
the impurity, Eqs. 2 and 3 reduce to A7iP*(r, T) = 
-|lm/(2r,r) • i-iyC, with f{r,T) = -inT/ smh{nrT) 
the Fourier transform of the Fermi function, and where 
G = 2er/(e2 F^) = sin(25) depends purely on 
the Fermi liquid phase shift 6. At low-temperatures 
T ^ l/r, one obtains asymptotically An^''{r,T) ~ 
— {—ly/r, while at higher temperatures, density oscil- 
lations are exponentially suppressed via An^'^{r,T) = 
4CT(— l)*" exp(— 27rr/<^T), with the thermal length-scale 
= l/T so defined, and appearing naturally. 
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Figure 3: Effect of temperature on the Anderson model {IQPV = 3, lO^e = —0.8, IQPU = 1). Panels (A) and (B): Comparison 
of impurity entropy with length-scales (solid lines) and ^lm (dashed) as they evolve with T. (C): |An(r, T)| vs r for 
T /Tk = 0.1, 0.5, 0.75 and 0.85. (D): £,k{T)/^k{0) vs T/Tk for the same system (circles), and scaling collapse of the divergence 
to common universal curve for V — 0.05, e = —0.1, U = 0.3 (crosses). Dashed line is a fit to the data. 



The situation of interest is of course the U > inter- 
acting AIM, where Gr,ai<-^,T) and hence T^{uj,T) now 
contain information about the Kondo effect. In this work 
we employ the numerical renormalization group (NRG) 
to obtain accurately Tc,{uj,T)}^ Provided T <^ Tk, sig- 
natures of Kondo physics and RG flow are naturally ex- 
pected in the Green function or LDOS for any site in 
real-space (see Eq. 3), and this has been recently con- 
firmed experimentally.^^ But the factorized form of Eq. 3 
implies that such measurements cannot reveal a Kondo 
length-scale, (there can be no universality in terms of 

The key point is that the real-space charge densities 
are obtained from a simple integral transformation of 
the frequency-resolved t-matrix by Eqs. 2 and 3. In- 
deed, far from the impurity this transformation reduces 
to a Fourier transformation, whose well-defined inverse 
guarantees the complete conservation of information for 
the process Tc{uj,T) o An{r,T). In consequence, the 
full RG structure of the underlying quantum impurity 
problem (as manifest in the t-matrix or LDOS) is wholly 
reproduced in the spatial variation of An(r, T). 

This fact is shown at zero-temperature in Fig. 2, where 
the energy-dependence of the spectrum —irp Im T„{Ld) 
is compared with the space-dependence of the densities 
An(r)/ AnP''{r), normalized by the pure potential scat- 
tering contribution at long distances (with C = sin(2(5) 
and 6 now the phase shift of the interacting problem). 
The representative systems plotted in Fig. 2 have similar 
energy scales Tlm ~ 0{U), but widely differing Kondo 



scales Tk- The densities shown in the lower panel exhibit 
clearly RG flow between FPs, now as a function of dis- 
tance: a non-universal crossover length-scale ^lm char- 
acterizes moment formation on flowing from FO to LM 
FPs, while (^k sets the universal length-scale for Kondo 
screening of this moment on flowing to the SC FP. 

Since the SC FP can be thought of as a free conduction 
band with one site removed, a 7r/2 phase shift results. -"^ As 
shown in Ref. 4 for the Kondo model, this is responsible 
for a sign-change in the Friedel oscillations. Indeed, this 
provides a natural explanation of the limiting universal 
values An{r)/AnP''{r) = -1 and -1-1 for r < (> ^lm) 
and r 3> ^iv; and the resulting minimum in | An(r)| allows 
direct identification of the Kondo length-scale. However, 
we point out that the FO FP can similarly be thought 
of as a free conduction band with an additional poten- 
tial scatterer. The resulting phase shift S also results in 
a minimum of |An(r)| at r = ^lm- Both length-scales 
of the Anderson model can be simply identified in this 
way, and we find ^k ^ 1/Tk and ^lm ^ 1/Tlm; consis- 
tent with expectations based on simple scaling grounds. 
In particular, grows as the impurity-host coupling 
decreases.^ Indeed, Kondo screening is totally absent in 
the limit of an uncoupled impurity y — > 0. Here local 
moment physics persists and ^k diverges. This implies 
directly that the spatial region r <C (^ Clm) is de- 
scribed by the LM FP, and cannot be as such a 'screening 
cloud'. 

One naturally expects the RG structure of the prob- 
lem to be manifest also in the temperature-dependence 
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of physical quantities. This is strikingly apparent in ther- 
modynamics such as the entropy, whose impurity con- 
tribution flows from 5inip(r) = ln(4) at the FO FP to 
ln(2) at the LM FP, and finally Sin,p = at the SC FP, 
symptomatic of complete Kondo screening; and as shown 
in Fig. 3(A). The corresponding effect of temperature on 
the real-space physics is demonstrated in panels (B)-(D). 
Two clear dips in | An(r, T)| are seen at each temperature 
in panel (C), allowing thereby immediate identification 
of £,lm{T) and £,k{T). As the temperature is increased 
through Tx, we find that (.k(T) moves out to larger dis- 
tances (while ^LM (T) remains unaffected) . On destroying 
the Kondo effect by increasing temperature, one would 
expect a 'screening cloud' per se to collapse; but under- 
standing instead as the length-scale for RG flow from 
LM to SC FPs, the divergence oi £^k{T) indicates simply 
the spatial persistence of local moment physics. Temper- 
ature cuts off RG flow to the SC FP, and so there can 
be no crossover Kondo length-scale to an SC region of 
space for Tr- < T < Tlm- The LM FP, describing an 
unscreened impurity, then pertains for r ^ ^lm- 

The characteristic divergence of S,k{T) exhibits uni- 
versal scaling in terms of T/Tk, as shown in panel (D). 
This behaviour is found to be described by £,k{T) ^ 
Cif (0)/[l - HT/Tk)] (dashed line), with b = C(l); and 
as such the diverging length-scale corresponds roughly 
to the vanishing energy scale (Tr- — T). The full evo- 
lution of ^liiT) and ^lmCT) is shown in Fig. 3(B). The 
free moment associated with the LM FP is destroyed by 
charge fluctuations at the FO FP when T ^ Tlm- ^lm 
thus diverges on warming through Tlm- The real-space 
behaviour in the three distinct temperature regimes is 
illustrated pictorially in Fig. 1. 
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Figure 4: T = spectrum vs inverse frequency (upper 
panel) and normalized density vs distance (lower panel) for 
the 2CK model with (Ji + J2) = 0.2. Symmetric case 
J\ ~ J2 shown as solid lines, while results for channel 1(2) 
shown as the dashed(dot-dashed) lines in the asymmetric 
case (Jl — J2) = 10~^. Common asymptotic behaviour for 
r < C2CK is Ani(r)/Anf (r) ~ c/ {r /^2Ck) - 1 while 
^ {r /^2Ck)~^'''^ behaviour is observed for r 2> ^2Ck («C C*)- 
In the asymmetric case, excellent agreement is seen with the 
exact result in the vicinity of the FL crossover, Eq. 6 (points) . 



III. MULTI-CHANNEL KONDO MODEL 

To emphasize the generality of the above results, we 
turn now to the A^-channel Kondo model, given by 

N 

K = Y. (^i-t + -h^ ■ §.) , (4) 

i=l 

describing a single spin-^ impurity S exchange-coupled 
to the conduction electron spin-density of channel i 
at the impurity. As before we consider a linear chain 
geometry, but with potential scattering^ in each chan- 
nel now included explicitly since the impurity is now 
strictly singly-occupied: Hl^^^ = Ea{A'4,o,<TC^,o,<T + 
YlT=a *(4,r,crCi,r-i-i,<T + H-c)}. Straightforward application 
of the equations of motion yield directly an analogue of 
Eq. 3 for each conduction channel i: 

AG,.,,,(w,T) = AGr(w) + [\X{Lo,r)Y T,,{uj,T), (5) 

where AG^''(aj) is the trivial potential scattering con- 
tribution in the Ji ^ {) case, X(aj,r) = Ar(w,r)/[1 — 
KG^ito)] contains the lattice information, and Tia-{u!,T) 
are the t-matrices of the multi-channel Kondo model 



(again obtained here using NRG). Friedel oscillations in 
each channel arise when K due to potential scat- 
tering [with C = sin(2(5) = -2Kt/{K^ + t^)]. But 
non-trivial many-body effects associated with the multi- 
channel Kondo model also appear in the real-space charge 
densities Ani{r,T) for Ji ^ 0, as now shown for the 
N = 2 two-channel Kondo (2CK) model.^O'^i 

In the channel-symmetric case Ji = J2 > 0, RG 
flow between two FPs determines the underlying physics. 
The LM FP describes the behaviour at high ener- 
gies/temperatures, where the impurity spin is essentially 
free. But intrinsic frustration arising from symmetric 
coupling to two conduction channels drives the system to 
the 2CK FP below T2CK, which describes a non- Fermi 
liquid (NFL) ground state characterized by 'overscreen- 
ing' of the impurity spin.^'''^^ This RG flow is of course 
clearly manifest in the behaviour of the t-matrix for each 
channel, Tia-{uj) (equivalent by symmetry), and so we ex- 
pect to observe the same flow in real-space. 

An analog to the single-channel 'Kondo screening 
cloud' has been suggested for the symmetric 2CK 
model: an 'overscreening' cloud of spatial extent 
£,2CK ^ l/72Cif surrounding a partially quenched im- 
purity. But since the ground state is a NFL, what lies 
beyond ^2CK? Within the present RG description, the 
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Figure 5: T = phase diagram for the 2CK model, showing 
LM, critical 2CK and FL regions of space. Direct crossover 
from LM to FL occurs for large asymmetry (diamond points), 
but a quantum critical region of space appears for smaller 
asymmetry. Two distinct scales are observed in this regime: 
^2CK (crosses) and 5* (circles). Solid line is fit to expected 
behaviour ^2Ck = 1/T2ck while dashed line is ^* = l/T*. 



impurity is surrounded instead by a LM cloud, with NFL 
behaviour described by the 2CK FP found in an over- 
screening cloud at distances r ^ ^20^- 

While the phase shift concept provides an intuitive ex- 
planation for physical behaviour in models with Fermi 
liquid ground states,^ we point out that this is inciden- 
tal. No such description in terms of phase shifts exists in 
the 2CK model, yet the full RG flow and universality of 
the problem is again recovered in real-space. 

The 2CK FP is in fact the critical point of Eq. 4.20.21 
For finite channel asymmetry ( Ji — J2) 7^ 0, a third FP 
enters: a Fermi liquid (FL) FP in which the impurity 
spin is completely screened by the more strongly-coupled 
channel. A new energy scale T* emerges, correspond- 
ing to RG flow from 2CK to FL FPs, and again shows 
up in the t-matrix for each channel. 2°"22 gj^^j 
pected that the corresponding length-scale ^* ^ l/T* 
characterizes real-space RG flow from a critical 2CK re- 
gion of space to an FL region. This is shown explicitly in 
Fig. 4 where we compare (as in Fig. 2) the spectrum vs 
inverse frequency (upper panel) to the normalized den- 
sities vs distance (lower panel). Two systems are plot- 
ted, with common K = 0.01 and (Ji + J2) — 0.2, but 
(Ji — J2) = for the symmetric case (solid lines), while 
(Ji — J2) = 10~^ for the asymmetric case (dashed lines 
for channel 1, dot-dashed lines for channel 2). The en- 
ergy scales T2CK and T* are indicated, as are the length 
scales ^2CK and ^* . 

Characteristic asymptotic behaviour in the vicinity of 
^2CK (given in the caption to Fig. 4, and shown as the 
dotted lines) can be extracted from Eqs. 5, 2 using the 
known asymptotics of the t-matrix for the symmetric 
2CK model. Indeed, the full density crossover from 2CK 
to FL FPs for r ^ £,2CK can be calculated using the ex- 
act t-matrix announced recently in Ref. 22. Our exact 



result at zero-temperature and small K follows as 



±Jf exp(r/r)V^Ko(r/r), (6) 



where Kg is the modifled Bessel function of the second 
kind, and the +(— ) sign is used for the density in channel 
1(2). The essentially perfect agreement between Eq. 6 
(points) and the full numerical solution (lines) is shown 
in the lower panel of Fig. 4. Interestingly, this result 
indicates that the density crossover approaching the FL 
FP of the 2CK model is closely related to the spatial 
crossover in the magnetization resulting from a boundary 
magnetic field acting in the Ising model. 22 

Finally, we plot the phase diagram for the 2CK model 
in Fig. 5. Three regions of space can be identified at zero 
temperature, corresponding to LM, critical 2CK and FL 
FPs. The crossover length-scales S,2CK ~ ^/{Ji ~ J2Y 
(dashed line) and £,2CK ^ exp[l/ max(/9Ji, PJ2)] (solid 
line) are anticipated from the well-known variation of 
the energy scales T2CK and T* in the asymmetric 2CK 
model, 2*''2i and agree excellently with explicit calcula- 
tion (points). Even in the 'standard' case where some 
degree of channel asymmetry is present in the model (in 
which case the ground state is a Fermi liquid), the non- 
Fermi liquid correlations at finite energies give rise to an 
intermediate region in space describing quantum critical 
behaviour. 



IV. CONCLUSION 

Quantum impurity problems are generically described 
in terms of an RG framework: as a function of energy or 
distance. Specifically, static charge density oscillations 
away from an impurity in real-space are simply related 
by integral transformation to the energy-resolved scat- 
tering t-matrix — a dynamic quantity, itself related to 
the impurity LDOS. As such, the well-known and rich 
behaviour associated with RG flow between fixed points 
is wholly reproduced in the densities. Indeed, this is not 
confined to the simplest 1-dimensional geometry consid- 
ered explicitly here. Although real-space behaviour does 
of course depend on the particular physical system under 
consideration, the underlying RG structure of the prob- 
lem must still appear. 

One implication is that the 'Kondo cloud' surrounding 
an impurity, typically defined as the region r < , ac- 
tually corresponds mainly to the LM FP. The screening 
process itself occurs on fiowing to the SC FP &i r ^ S^k- 
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